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We calculate the relaxational dynamical critical behavior of systems of 0(n||) © 0{n±) symmetry 
by renormalization group method within the minimal subtraction scheme in two loop order. The 
three different bicritical static universality classes previously found for such systems correspond to 
three different dynamical universality classes within the static borderlines. The Heisenberg and 
the biconical fixed point lead to strong dynamic scaling whereas in the region of stability of the 
decoupled fixed point weak dynamic scaling holds. Due to the neighborhood of the stability border 
between the strong and the weak scaling dynamic fixed point corresponding to the static biconical 
and the decoupled fixed point a very small dynamic transient exponent, of lo^ = 0.0044, is present 
in the dynamics for the physically important case ny = 1 and n± = 2 in d = 3. 

PACS numbers: 05.50.+q, 64.60.Ht, 64.60.-1 



I. INTRODUCTION 



namical models is given. 



The phase diagram of systems with 0(7i||) © 0{n±) 
symmetry contains several phases meeting in a multi- 
critical point. In Ref.^ (henceforce called part I) it was 
shown that the static critical properties can be quan- 
titatively analyzed from field theoretic functions in two 
loop order if one uses resummation. As an example we 
have in mind an antiferromagnet in an external magnetic 
field (with ri|| = 1 and nj_ — 2), although other physical 
examples with different values of order parameter (OP) 
components may be considered. 

In order to get more insight in the dynamical critical 
properties near such a multicritical point we reconsider 
the simplest dynamical model possible for 0(n||)©0(n_L) 
symmetric systems. In such a dynamical model one as- 
sumes relaxational behavior for the two OPs 0|| and 
This model has been briefly studied^ on the basis of the 
static one loop results^. Meanwhileii^i^ it has been shown 
that the one loop results^s^ are considerably changed in 
higher loop order concerning the regions of different static 
multicriticality in the space of OP components n|| and 
n±. For integer order parameter components only a sys- 
tem n|| = 1 and nj_ — 2 belongs to the universality class 
characterized by the biconical fixed point (FP) indicating 
tetracritical behavior - if the physical system lies in the 
attraction region of the FPX 

The paper is organized as follows: In chapter |TT] we 
define the dynamical model, then in chapter IIIII the dy- 
namical field theoretic functions are introduced and the 
results in two loop order are presented. From these re- 
sults the FP and dynamical exponents are calculated in 
chapter IIVI and the stability of the FP is considered in 
chapter |Vl Due to the small dynamic transient exponent 
found, the effective - nonasymptotic - dynamical behav- 
ior is studied in detail in chapter IVll Finally conclusions 
and an outlook to subsequent research of extended dy- 



II. DYNAMICAL MODEL 

The results obtained in part I for the statics of systems 
with 0(fi||) ® 0{n±) symmetry are applied to the criti- 
cal dynamics if the system dynamics is described by two 
relaxational equations for the OP components (/)_lo smd 
(/)||o in the two subspaces. Correspondingly two kinetic 
coefficients r± and Tn have to be introduced. The model 
A type equations are 
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The stochastic forces 9^j^ and 9^^^ fulfill Einstein relations 

{9^^ix,t) 9^^^ix',t'))=2t^S{x-x')S{t-t')d"'' , (3) 
(0;, (x, t) (x', i')) =2f ||<5(x - x')Sit - t')5^= , (4) 

with indices a, (3 ~ 1, • ■ • , n±_ and i,j = 1, . . . , 7i|| cor- 
responding to the two subspaces. The static functional 
H-Bi is defined as 
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The properties and renormalization of the static ver- 
tex functions following from Ti-Bi have already been pre- 
sented in part I. There in resummed two-loop approxima- 
tion it has been shown that within a small region in the 
space of the spatial dimension and the OP components 
the biconical FP is stable (e.g. for n\\ = 1 and = 2 at 
d = 3). 



III. RENORMALIZATION, FIELD THEORETIC 
FUNCTIONS 

From the dynamic equations ([T]) and ^ a functional 
may be derived which allows the calculation of dynamic 
vertex functions in perturbation theory (for an overview 
seefi"). Within this dynamic functional additional auxil- 
iary densities 4>^_q and 0||q are introduced^ . Recently it 
has been shown^P that the dynamic two point functions 
have a general structure, which is in the current model 



(6) 



+fj°)(e^,^ll,fc)f|l , (7) 



where f 
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_L± (^-LjC||jfc) ^'^'^ ^iiii' (C-i-:'f|h^) are the static 
two point vertex functions discussed in part I. The func- 
tions (l^^^^{^^,i\\,k,u) and ^^^^^^^ (^±,C||,fc,w) have to 
be determined within dynamic perturbation expansion. 
All functions in ^ and (O depend besides the correla- 
tion functions , , the wave vector modulus k, and the 
frequency also on the static couplings Ux and 
The functions Cl^ ^ (with at =_L, ||) additionally de- 
pend on the two kinetic coefficients fx and fji . As we will 
see below, the genuine representation ([6]), ([7| that allows 
to single out contributions from merely static vertex func- 
tions into dynamic ones essentially simplifies cumbersome 
calculations and enables one to effectively proceed with 
calculation of the dynamic RG perturbative expansions. 



A. Renormalization of tlie dynamic parameters 

The renormalization of the static quantities appearing 
in ([5]) has been presented in part I in detail and explic- 
itly performed in the minimal subtraction RG schemeii 
directly at d = 3 to the two-loop order. The resulting 
renormalization factors and field theoretic functions {C,- 
and /3- functions) remain valid also in dynamics. Addi- 
tional renormalizations are necessary for dynamic quan- 
tities. Within the current dynamic model only the auxil- 
iary densities and the kinetic coefRcients have to be renor- 
malized. 



The renormalized counterparts of the auxiliary densi- 
ties are defined as 



= z 



1/2; 



'110 - n ■ (8) 

The renormalized kinetic coefficients are introduced as 

fx = Zr^rx, f||=Zr„r||. (9) 

Relation ([5]) and the renormalization of the OP densities 
0j_o and <^||o introduced in part I imply for the dynamic 
vertex functions the renormalization 



- <t>\\(, 
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(10) 
(11) 



From the above relations and the structure of the dy- 
namic two point vertex functions presented in ([6]) and 
([7]) follows that the renormalization factors of the kinetic 
coefficients P^ and Pn in case of the absence of mode cou- 
plings are determined by the corresponding renormaliza- 
tion factors of the auxiliary densities. This leads to the 
relations 



Zr. = Z'Z 



l/2„-^l/2 



4'±- tj}± 



Zv.. = ZJ Z, 



1/2 ^7-1/2 
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The static renormalizaton factors Z^^ and Z^^^ have been 
introduced in Eq. (4) of part L 

B. Dynamic 13- and (^-functions in two loop order 

Quite analogous to statics in part I we will use the 
uniform definition 



Ca.(M,P±,P||) = 



din K 



(13) 



for the (^-functions also in dynamics, where is now 
a placeholder for any auxiliary density or kinetic coeffi- 
cient, K is the scaling parameter, and {u} = {u^, Ux , u\\ } 
is the set of static couplings. From perturbation expan- 
sion the resulting two loop expressions for the ^-functions 
of the kinetic coefficients P 1 and Pii read 
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The important dynamic parameter is the time scale ratio 



(16) 
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between the two kinetic coefficients T± and Fn, which 
have been aheady introduced in and p^ . From 
the above definition of the time scale ratio and the def- 
inition of the C-functions in (|13p the /3-function of v is 
determined by 



f3v 



dv 
dn 



(17) 



where the derivative is taken at fixed unrenormalized 
quantities. Inserting (|14p and (fT5|) into (fT7|) the two loop 
expression of the /3-function of v reads^'^^ 
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The /3-function changes its sign under interchanging the 
parallel and perpendicular components and replacing the 
time scale ratio v by 1/v. 

In the nonasymptotic region where a non universal ef- 
fective critical behavior may be observed the values of the 
static couplings and the time scale ratio are described by 
the flow equations. For v it reads 



= /3„(u||(/),Ux(0>"x(0) : 



(19) 



whereas for the static couplings Eq. (36) of part I with 
the Borel resummed static /3-functions are used. The 
asymptotics is reached in the limit I —^ starting in the 
background at / = 1 from non universal initial values of 
the time scale ratio and couplings. 



FIXED POINTS AND DYNAMICAL 
CRITICAL EXPONENTS 



IV. 



As usuals the two C- functions Eqs. p4|) and (|T5l) define 
two dynamical critical exponents that govern the power 
law increase of the autocorrelation time for the OPs chn 



and (/)_L, correspondingly 

z± = 2 + and 



2 + C?„ 



(20) 



where the stable FP values of the static and dynamic 
parameters have been inserted into the ^-functions, this 
means Cp = Cr„. ({w"^}, f*) . At the strong scaling FP 
there is only one dynamic time scale and the two ex- 
ponents are equal whereas at the weak scaling FP they 
are different and define for each component, parallel and 
perpendicular, the time scale. 

Depending on the FP value of the time scale ratio v 
one may obtain strong (v* ^ 0, oo) or weak {v* = 0, oo) 
dynamic scaling. The dynamical FPs are calculated (see 



also Eq. (12) in Ref.^) from setting the /3-function (fT8)) 
equal to zero. Inserting the stable static FP values (see 
Table I in part I) into Eq. (|18p one then may calculate 
a dynamical 'phase diagram' in the 7i||-rtj^-plane quite 
similar to the static 'phase diagram' Fig. 1 in part I. 
Let us note here, that one can make use of two different 
ways to analyze perturbative expansions within the mini- 
mal subtraction RG scheme. The first one is the familiar 
e-expansion, when the FP coordinates and asymptotic 
critical exponents are obtained as series in £ = 4 — c? 
and then evaluated at the dimension of interest (e.g. for 
d = 3). The second one relies on treatment of the expan- 
sions in renormalized couplings directly at fixed dimen- 
sion d = 3}^ Enhanced by resummation such a scheme 
allows to treat, besides the asymptotic quantities, the 
non-universal effective exponents. The latter method has 
been applied in part I to perform a comprehensive analy- 
sis of non-universal static behavior. Below we will make 
use of the static results obtained there to proceed with 
the analysis of (asymptotic and effective) dynamical crit- 
ical behavior. 

To summarize an outcome of the static FP stability 



analysis 



1.3.4.5.6 



let us recall that, depending on the nu, 



values the critical behavior is governed by one of 
the three non-trivial FPs: (i) isotropic Heisenberg FP 
H(n_L -f ny) with My = = u* — u*; (ii) decoupling FP 
V with u| ^ 0, 7^ 0, = 0; (iii) biconical FP B with 
U|l 7^ 0, ^ 0, 7^ 0. Below, we will analyze pecu- 
liarities of the dynamical critical behavior in the above 
universality classes. 



A. Dynamics at the isotropic Heisenberg fixed 
point 

At the isotropic Heisenberg FP H{n± +n||) the fourth 
order static couplings are equal, u| = = u* = u*. 
In consequence the static couplings drop out in the FP 
equation for v. Assuming a nonzero finite value of v at 
the FP, the equation for v*, Pviv*) = 0, reads 



(nil + 2) - (n^ + 2) (6 In 
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One immediately sees that for general and n± a zero 
can only be found if the arguments of the logarithms are 
equal to 4/3, which leads to the FP value 



1 



(22) 



This result has to be fulfilled also in higher loop order 
and therefore the result is exact. Due to Eq. (|^^ the 
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(■-functions ([14]) and (fTSl) become equal at the FP 

C?, =C?„ =(" + 2K(61n^-l) (23) 

with n = ny + n_L. This means strong dynamic scaling 
with the dynamical critical exponent 

z = 2 + cri (24) 

of the 0(n)-symmetric model A universality class, where 
c = 6 In I — 1 in two loop order,— and 77 is the anomalous 
dimension of the 0(n)-synimetric model. 



B. Dynamics at the decoupling fixed point 

At the decoupling FP V static critical behavior does 
not fulfill scaling due to the existence of two different 
correlation lengths, each for one of the decoupled parts of 
the system (the parallel and perpendicular components 
of the OP). The FP value of the static coupling Ux is 
equal to zero. Therefore the flow equation p9)) at the 
static FP reduces to 

1^ = ^ ([(nil + 2)uf [n^ + 2)<^] (6 In ^ - 1)) (25) 

leading either to a flow reaching v* — or \/v* — 
depending wether 

n|| < nj_ where Cr^'^'' Cr*^"^' or 

nil > where Ct\^'^ > C*^^^ , (26) 

with Cfj"^^ ft^s (-function of model A for the correspond- 
ing subsystem with n|j_L components. Both cases mean 
that weak scaling holds at the decoupling FP. Indeed 
inserting the values of the decoupling FP into the (- 
functions and gives two dynamical exponents 
z. One for the dynamics of the parallel and another for 
the perpendicular components of the OP. Both exponents 
correspond to the model A universality class 

Z|| = 2 + c?7|| and = 2 + cry^ . (27) 

In the special case when n|| = nj_ the exponents Z|| and 
z± become equal and the FP values of the time scale ratio 
are determined by the initial values of the static flow. 



C. Dynamics at the biconical fixed point 

At the biconical FP B the FP values of the three fourth 
order couplings are different and the solution for the FP 
value of V becomes nontrivial and dependent on the num- 
ber of components n|| and n±. However the only relevant 
case where this FP is stable in c? = 3 is (ny, n±) = (1, 2) 



or (n|| , n±) = (2, 1) j2, Due to symmetry the two solutions 
are related 

^;*(n||,nj_) = — i ^. (28) 

The numerical solution found using the static one loop 
order FP values for the static couplings reads^ 

w*(l, 2) = = 1.0241. (29) 

Thus one finds strong scaling with a new biconical dy- 
namical exponent. Inserting the FP values into the (- 
functions HH) and the dynamical critical exponent 
is given by 

= 2.0149. (30) 

As has been shown in part I the biconical FP becomes 
stable in two loop order within a small region around the 
OP-component values nj_ — 2 for n\\ =1. However one 
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FIG. 1: Dependence of FP values of the timescale ratio Ijv* 
and the dynamic transient exponent ld^ on tix for n\\ = 1 in 
the region of stability of the Heisenberg 7i(3) and biconical B 
FP in d = 3. The dotted vertical line indicates the stability 
border between 7i(3) and B. At the stability border to the 
decoupled FP both Ijv* and u)^ go to zero (at nx ~ 2.18). 

has to apply resummation techniques to the two loop 
functions in order to get real FP values for the static 
couplings. Using these two loop order resummed values 
the dependence of the FP value of the time scale ratio v 
within this region is shown in Fig. [TJ Note that we do not 
resume the expression for /3„ , Eq. (fTS]) , itself. The biconi- 
cal FP reaches the value of the Heisenberg FP, — 1, 
at the stability border line and the decoupled FP value, 
?;* = or w* = 00 (depending on wether n|| is larger or 
smaller than n±) at the corresponding stability border 
line. Inserting for ny = 1 and n± = 2 the resummed FP 
values for static fourth order couplings into Eq. one 
obtains 

V* = v'^ = 1.0555. (31) 
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The corresponding dynamical critical exponent reads 



= 2.052. 



(32) 



This has to be compared with the predicted value for the 



Heisenberg FP- 



2.015, which was found to be 



stable in the e-expansion in one loop order . 

V. DYNAMIC TRANSIENT EXPONENTS 

The static stability boundaries are also dynamic sta- 
bility boundaries and therefore in the case ny = 1 and 
nj_ = 2 a small dynamic transient exponent is expected. 
Its value is given by 



dv 



(33) 
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It will be further numerically evaluated by inserting the 
Borel resummed values for static couplings. As will be 
shown below this exponent goes to zero only when the 
dynamical FP changes from the strong dynamic scaling 
to the weak dynamic scaling FP. This is the case when in 
addition to the change of the stability of the static FP, 
the stability of the dynamic FP is changed. 

The instability of the weak scaling FP w* = or 1/t;* = 
is defined by a negative dynamic transient exponent 




(61n--l) 



nil u 



*2 



2 lim ln(2i;) 



(34) 



This shows that in any case where u*^ is different from 
zero the weak scaling FP is never stable. Thus only at 
the decoupled fixed point the dynamic transient exponent 
might be positive, if ny > n±, see Eq. (|26p . The transient 
exponent then reads 



OJ.y 



18 



[V\\ - V±) 



(35) 



At the Heisenberg FP the dynamic transient exponent 
reduces to 



n 



n(w^("))2 , 4 

In — 

18 3 



(36) 



Thus the dynamic transient exponent at the stability bor- 
derline to the biconical FP is finite and continous (see 
Fig. [T]). In the region of stability of the biconical FP 
the dynamic transient exponent is given by the expres- 
evaluated with the appropriate FP values for Ux 




weak scaling 
z =2+cri 

par 'par 

Z =2+Cri 
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FIG. 2: Regions of different static bicritical behavior, which 
are defined by the stable FP in the ny — n^-plane (e = 4 — d = 
1) separated by the static stability border lines, (from left to 
right: Heisenberg FP, biconical FP and decoupled FP). The 
solid line is also the border line between dynamic strong and 
weak scaling. The dots indicate small integer values of the 
component numbers. 



and V. At the stability border lined to the decoupling FP 
both FP values go to zero. Thus also the dynamic tran- 
sient exponent goes to zero indicating the change from 
the stability of the strong scaling dynamic FP to the weak 
scaling FP. Inserting the FP value for the biconical FP 
leads to a dynamic transient exponent roughly one order 
smaller than at the Heisenberg FP due to the smaller FP 
value of the static coupling Ux- Inserting the FP values 
for the biconical FP into (l33l) one obtains 



= 0.0044 . 



(37) 



Thus in addition to the already small transient from stat- 
ics an even smaller transient in dynamics appears. This 
leads to a slow approach of the FP values in the flow 
equations. 

The resulting 'phase diagram' concerning the dynam- 
ical universality classes is shown in Fig. [21 A strong 
dynamic scaling part at small values for the OP compo- 
nents is separated by a stability border line (solid curve) 
at which the dynamical transient uj^ goes to zero. This 
border line lies very near the dot representing the model 
describing the critical behavior of a three dimensional 
Heisenberg antiferromagnet in a magnetic field (ny = 1, 
nj_ = 2). In consequence the transient from the back- 
ground to the asymptotic behavior might be very slow. 



VI. FLOW EQUATIONS AND EFFECTIVE 
EXPONENTS 



The asymptotic dynamic exponents may be reached 
only in very small region around the FP where the devi- 
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FIG. 3: Effective dynamic exponent for the static flows 3, 5 
and 6 (see Fig. 1 in part I). The timescale ratio v is set to its 
biconical FP Eq. (|3T| . 



FP both exponents reach the same value since the strong 
scaling dynamic FP is stable. In order to show the effect 
of the small static transient exponent we fix the value of 
the time scale ratio to its biconical FP value. The ini- 
tial values for the static couplings are chosen to be the 
same as in part I for the flows in Fig. 1 numbered from 
1 to 6. Although the static FP value is not reached (see 
Fig. 1 in part I) the numerical differences in the effective 
dynamical exponents are small. The difference between 
the parallel and perpendicular effective dynamical expo- 
nent for curves number 2 and 3 in the background region 
(larger /) result from that part of the static flow where 
Ux and either M|| or u± are almost zero. 

In order to show the effect of the smaller dynamical 
transient exponent ([37]) we fix the static couplings to 
their biconical FP values and start the flow for the time 
scale ratio v at three different initial values correspond- 
ing to situation where the parallel relaxation coefficient 
is smaller, equal or larger to one. 



ation from the FP values in the model parameters have 
died out. Due to the small transient exponents (either 
static and/or dynamic) in the physical accessible region 
the critical behavior may be an effective one described by 
effective exponents calculated with the parameters differ- 
ent from their FP values at I — and obtained from the 
flow equations at finite values of I. The effective expo- 
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FIG. 4: Effective dynamic exponent for the static flows 1, 2 
and 4 (see Fig. 1 in part I) . The timescale ratio v is set to its 
biconical FP Eq. (|3T|l . 

nents are defined as 

ZeffAl) = 2 + Cr, (ui(/), u|| (0, (0, ^^(0) , (38) 

ze//,||(0 - 2 + Cr, Ml), u\i{l),ux{l),v{l)). (39) 

In Figs. [3]and|3]we show the effective dynamic exponents 
for the parallel and perpendicular components of the OP. 
In the asymptotics when reaching the stable biconical 
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FIG. 5: Flows of the timescale ratio v(l) at the biconical 
static FP B for different initial conditions. v{l = 1) = 0.1, 
v{l = 1) = 1 and v{l = 1) = 5. Note the scale of the flow 
parameter compared to Figs. I3I4I6I 

As one can see from the Fig. [S] indeed for initial values 
of v{0) far from its FP value, the timescale ratio almost 
never attains its FP value (for v{l) = 0.1 it reaches the 
asymptotics at In^ ~ —10^)! However the effective expo- 
nents Zeff are not so far from their FP values in conse- 
quence of the general dependence on the timescale ratio 
(see Fig. [6|). However one might define a dynamic am- 
plitude ratio from the relaxation rates Fn/Fx. This ratio 
then would in leading order behave like v. 

Starting the flows of 3 and 5 with v different from its 
FP value (see Fig. [7]) leads to nonmonotonic behavior of 
the effective dynamic exponents (see Fig. At the 

first sight the behavior of the parallel and perpendicu- 
lar effective dynamical exponents looks strange since the 
nonmonotonic behavior is seen in Fig. [5] for small ini- 
tial values v(0) in the perpendicular exponent, whereas 
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FIG. 6: Effective dynamical exponents z^jj,\\ (solid curves) 
and Zsff,± (dashed curves) for the flows shown in Fig. [S] 




FIG. 8: Effective dynamic exponent z^ff^w (solid curves) and 
Zeff,± (dashed curves) for the static initial conditions of flow 
number 5 and dynamic initial conditions as indicated. For the 
different behavior of the parallel and perpendicular exponent 
see the text. 



depending on v and the static coupling u^- For flows 
where Ux is very small no difference is seen according to 
V. So no difference is seen in the corresponding flow for 
the static flow number 3. 

However for Ux near the bicritical FP, the difference 
between and z± calculated at the same values of the 
static couplings but at different values of the timescale 
ratio sometimes can be positive (i.e. z\\ > z^) and some- 
times it can be negative (^n < zj^), depending on par- 
ticular values of v. Numerically estimates indeed recover 
the differences shown in Fig. [5] 



FIG. 7: Flows of the timescale ratio v{l) to the biconical 
static FP B for different initial conditions. 3: v{l = 1) = 0.1, 
5: v{l = 1) = 5. A quite different behavior is observed. 



VII. CONCLUSION AND OUTLOOK 



for large values v{Q) in the parallel exponent. To explain 
such an unexpected behavior, one may look at the dif- 
ference of the two effective exponents (see Eqs ([T4l and 

USD), 

Az = Ze^/j|(u||,Mx,U±,w) - Ze//,±(W||,Ux,U_L,l') (40) 

and estimate it at fixed U|| , Ux j '^.l but for different values 
of V. The difference can be written as: 



Az = c + 5z{v). 



(41) 



In Eq. (|4ip . c depends on the static coupling only and 

5z{v) for n±_ = 2, rt|| = 1 reads: 



^ ' 36 I V 1 + 2?; 
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(42) 



We have reconsidered the relaxational dynamics at the 
multicritical dynamical FPs in 0{n\\)(BO{n±) symmetric 
systems. According to the static two loop order results 
the biconical FP is the stable FP at the interesting case 
n|| = 1 and n± — 2 for which a new dynamic FP with 
strong dynamic scaling has been found. 

The critical dynamics of such a system has to take into 
account additional properties, namely if the densities of 
conserved quantities couple statically to the OP and/or 
if mode coupling term are present. Both extensions of 
the dynamical equations have to be considered lacking a 
complete two loop order calculation^. The model C like 
extension of this model will be presented in a third part 
of this series. 
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